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Ionization  Gauge  Measurements  of  Atmospheric 
Density  From  a Low  Altitude  Satellite 


l.  INTRODUCTION 

One  of  the  experiments  flown  onboard  an  Air  Force  satellite  launched  in  Oct- 
ober 1974  was  a cold  cathode  ionization  gauge.  Measurements  were  acquired  over 
a 7-month  period,  the  life  of  the  satellite,  into  May  1975.  Atmospheric  densities 
deduced  from  the  measurements  performed  in  the  period  from  November  1974  to 
January  197  5 have  been  processed  and  are  given  in  Appendices  A and  B.  Mea- 
surements performed  after  January  1975  are  currently  being  processed  and  will  be 
presented  in  a future  report.  The  density  results  are  for  an  altitude  region  near 
perigee  covering  a symmetric  altitude  range  centered  around  the  perigee  and 
extending  in  height  up  to  approximately  300  km  (see  Figure  1). 

Depending  upon  the  location  and  altitude  of  perigee,  a recorded  perigee  pass 
could  vary  in  time  from  15  to  20  min  and  cover  a latitude  span  of  approximately  70 
deg.  Flight  operations  and  the  data  acquisition  plan  are  more  fully  described  in 
Section  2 of  this  report. 

Other  sections  include  a description  and  discussion  of  the  satellite,  the  experi 
ment,  theory  of  measurement,  and  the  processing  and  reduction  technique  used  to 
obtain  the  results.  Density  results  are  given  in  Appendices  A and  B.  The  re- 
sults are  presented,  however,  without  any  attempt  to  analyze  or  interpret  them  in 
terms  of  atmospheric  structure,  as  that  is  planned  for  a future  date. 

(Received  for  publication  24  May  1976) 
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Figure  1.  Orbital  Segment  Depicting 
Perigee  Pass.  Data  acquisition 
occurs  during  the  10-min  periods 
before  and  after  perigee 


2.  SATELLITE  AND  FLIGHT  OPERATIONS 

The  primary  objective  of  the  Air  Force  satellite  mission  was  to  investigate 
and  perform  measurements  of  thermospheric  properties  that  influence  upper 
atmospheric  structure.  To  accomplish  this  objective,  the  spacecraft  was  instru- 
mented with  a complementary  payload  consisting  of  density,  composition,  and 
energy-source  experiments. 

The  spacecraft  payload  was  placed  into  an  elliptical  polar  orbit  inclined  at 
96.  98  deg,  with  perigee  and  apogee  heights  of  158.  6 and  3761.  8 km  respectively. 
Initial  perigee  location  was  at  44.  9 deg  north  latitude,  and  it  precessed  northward 
at  a rate  of  2.  22  deg  per  day.  Following  injection  of  the  satellite  into  orbit,  the 
satellite  was  maneuvered  into  its  operational  attitude,  which  was  spin  stabilized 
and  oriented  to  have  the  spin  axis  directed  within  5 deg  of  the  normal  to  the  orbit 
plane. 

The  maneuvering  and  maintenance  of  the  satellite's  attitude  were  accomplished 
by  means  of  magnetic  torquing  coils  that  interacted  with  the  earth's  magnetic  field 
to  provide  the  forces  needed  to  control  or  maintain  the  desired  orientation.  The 
magnetic  torquing  system  functioned  smoothly  throughout  the  satellite's  lifetime. 

The  adjunct  to  attitude  control  is  attitude  determination,  and  this  information 
was  provided  by  a combination  of  earth  and  sun  sensors  whose  outputs  served  as  a 
primary  source  of  reference  data.  Backup  capabilities  were  provided  by  a three- 
axis  magnetometer.  This  backup  capability  was  utilized  very  early  in  the  flight  as 
the  earth  sensor  responded  erratically  when  the  vehicle  was  sun-lit.  The  opera- 
tional attitude  was  then  determined  from  a sun  sensor-magnetometer  combination 
during  daylight  hours,  and  from  an  earth  sensor-magnetometer  combination  when 
the  satellite  was  occulted. 

Figure  2,  a photograph  of  the  satellite,  shows  the  spacecraft's  magnetometer 
boom,  some  solar  panels,  and  some  payload  sensors  mounted  in  their  flight  con- 
figuration. The  earth  sensor  can  be  seen  located  about  midpoint  on  the  surface. 
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Figure  2.  Spacecraft  Shown  With  Solar  Panels  and  Thermal  Blankets  Installed. 

On  the  forward  surface  some  experimental  sensors  can  be  seen  in  their  flight 
configuration 

at  the  right  of  the  figure,  embedded  in  the  solar  panels.  The  solar  panels  were 
the  principle  power  source,  operating  satellite  systems  directly  during  daylight 
hours.  For  night  operation,  systems  were  powered  by  batteries  that  had  been 
charged  from  the  solar  cells.  Normal  on-orbit  support  consisted  of  payload 
operation  during  alternate  orbits.  Where  a payload  operation  consisted  of  recorded 
measurements  over  a perigee  pass  (which  is  an  orbital  segment  centered  at  peri- 
gee), the  pass  was  automatically  cycled  to  occur  at  a predetermined  height  before 
perigee  and  to  terminate  at  the  same  predetermined  height  after  perigee.  Once 
per  week,  in  addition  to  perigee  pass  coverage,  a full  orbit  of  payload  operation 
was  recorded. 

In  the  event  of  an  atmospheric  disturbance— for  instance,  magnetospheric 
storms  or  strong  solar  flare  activity,  a plan  for  additional  support  was  formulated. 
Criteria  had  been  established  that  involved  setting  up  alert  indices  with  threshold 
values  that  when  exceeded  would  identify  the  onset  of  the  activity  event.  With  the 
triggering  of  the  threshold  value,  additional  support  coverage  was  initiated  and 
continued  as  long  as  the  event  was  in  progress.  Other  criteria  had  been  established 
to  identify  when  the  event  had  subsided;  once  this  occurred,  special  coverage  would 
cease  and  support  would  revert  to  normal  orbital  operations.  In  actual  practice. 
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the  additional  support  plan  meant  that  a maximum  number  of  perigee  passes  con- 
sistent with  the  spacecraft's  power  and  thermal  system's  constraints  and  limita- 
tions would  be  obtained.  Both  plans,  for  normal  and  special  coverage,  functioned 
without  a major  disruption  throughout  the  satellite's  lifetime. 


3.  INSTRUMENTATION 

Figure  3 shows  the  ionization  gauge  experiment  that  was  flown  on  the  Air  Force 
satellite.  The  unit  shown  on  the  right  of  the  figure  houses  the  electronic  circuitry 
required  for  the  operation  of  the  gauge  sensor.  Included  in  this  unit  are  high  and 
low  voltage  power  supplies,  monitor,  calibration  and  conditioning  circuits,  and  a 
linear  automatic- range-switching  amplifier.  The  unit  with  the  reflective  front 
surface  is  the  sensor  unit;  it  is  shown  with  the  gauge  uncapping  mechanism  in  place, 
as  it  is  before  activation.  Figure  4 shows  the  ionization  gauge  sensor  and  its  mag- 
net after  calibration  but  before  installation  into  a flight  package.  The  ionization 


to  thrift  IOrLZatiT  <t,a'ige  FJight  Unit:  Electronics  are  Packaged  in  the  Unit  Shown 
to  the  Left.  The  polished  surface  area  on  the  sensor  unit  is  the  surface  exposed  to 
the  ambient  atmosphere  r 
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Figure  4.  Cold  Cathode  Gauge 


gauge  is  a cold  cathode  magnetron  type.  The  connection  with  an  attached  demount- 
able flange  located  to  the  left  of  the  gauge  body  is  fcr  calibration  purposes.  Follow- 
ing calibration,  the  gauge  is  evacuated  and  sealed  preparatory  to  being  assembled 
into  a flight  unit.  The  sealing  off  process  is  performed  by  "pinching  off"  the  copper 
calibration  stem,  which  can  be  seen  to  the  left  of  the  flange  attachment.  Extending 
from  the  right-hand  side  of  the  gauge  body  is  the  sampling  tabulation.  This  is  the 
connection  that  is  opened  to  the  atmosphere  once  the  gauge  has  been  placed  in  orbit. 
The  gauge  opening  or  uncapping  mechanism  fits  about  the  sampling  tube  and  is 
firmly  attached  to  the  tube  by  means  of  a lever  arm  that  is  fastened  onto  the  tube 
end  cap.  Upon  activation  of  the  lever  arm,  a torque  is  created  that  ruptures  the 
sampling  tube  and  exposes  the  gauge  to  the  environment.  The  break-off  interface 
for  the  sampling  tube  is  defined  by  the  groove  located  just  to  the  rear  of  the  end 
cap  (Figure  4).  At  the  midpoint  of  this  groove  the  wall  thickness  has  been  reduced 
from  0.  125  in.  to  0.  005  in. 

The  gauge  is  opened  to  the  ambient  atmosphere  by  means  of  explosive  bellows 
that  are  enabled  and  fired  by  ground  command  some  24  hr  after  orbit  has  been  ob- 
tained. The  gauge  end  cap  together  with  the  break-off  mechanism  are  ejected. 
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leaving  a 2 rr  steradian  viewfield  at  the  sampling  aperture.  The  delay  in  opening 
and  operating  the  gauge  allows  time  for  the  atmospheric  gases,  brought  aloft  by 
the  booster  and  spacecraft,  to  vent  into  space  - a precaution  taken  to  guard  against 
the  occurrence  of  corona  inside  the  payload. 

Normal  data  acquisition,  as  described  in  Section  1,  begins  after  evaluation 
and  confirmation  of  experiment  operation.  The  gauge  sensor  located  in  the  spin 
plane  on  the  spin  equator  rotates  with  the  satellite  so  that  the  sampling  aperture 
alternately  faces  into  and  away  from  ram.  As  a consequence  of  this,  the  gauge 
internal  pressure  is  modulated,  increasing  to  a maximum  at  ram  and  subsequently 
decreasing  to  a minimum  at  wake.  Ram  is  that  portion  of  the  spin  cycle  where  the 
angle  between  the  velocity  vector  and  the  sampling  tube  axis  is  less  than  90  deg, 
with  maximum  ram  occurring  at  a zero  degree  attack  angle.  Within  the  gauge,  the 
internal  gas  is  ionized  and  a plasma  forms  that  is  confined  to  the  discharge  region. 

Positive  ions  extracted  from  the  plasma  produce  a current  output  that  varies 
in  magnitude  in  proportion  to  gauge  pressure.  The  gauge  internal  pressure,  in 
turn,  is  a function  of  ambient  density  and  attack  angle,  defined  as  the  angle  between 
the  velocity  vector  and  the  axis  of  the  gauge  sampling  tube.  A typical  modulated 
gauge  response  is  shown  in  Figure  5.  The  angular  resolution  of  the  gauge  response 
is  determined  by  the  sampling  rate.  That  rate  for  this  flight  was  16  samples  per 
sec,  which  was  approximately  1 measurement  every  2 deg  of  rotation. 


4.  THEORY 
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A number  of  investigators  ’ ’ have  developed  the  theory  relating  ambient 
density  to  ionization  gauge  measurements  under  the  conditions  of  satellite  flight. 
The  expression  relating  ambient  pressure  to  gauge  pressure  is: 


P 

g 


R(S,  D,  a) 


<1' 


where 


1.  Pond,  H.  L.  (1962)  The  effect  of  entrance  velocity  on  the  flow  of  a rarefied  gas 

through  a tube,  J.  Aerospace  Sci.  , £9:917-920. 

2.  Ivanoskii,  A.  I.  (1964)  Interaction  between  an  instrument  cavity  and  a flow  of 

rarefied  gas,  Trudy,  5^:49-96. 

3.  Hughes,  P.  C.  (1965)  Theory  for  the  Free  Molecular  Impact  Probe  at  an  Arbi- 

trary Angle  of  Attack,  University  of  Toronto  UTIA.S  Report  No.  f63. 
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GAUGE  PRESSURE  (Torr) 


Figure  5.  Gauge  Response 
Curves  Showing  Variation  of 
Gauge  Internal  Pressure 
Over  Two  Satellite  Spin 
Cycles.  Maximum  ram 
points  are  at  altitudes  of 
215  km  and  208  km 


R(S,  D,  Of  1 

S 

u 

v 

m 


gauge  internal  pressure  and  ambient  pressure,  respectively 

absolute  temperature  of  the  gauge  and  ambient  gas,  respectively 

gauge  sampling  tube  factor,  a function  of  transmission 
probabilities 

speed  ratio  u/v 

m 

satellite  velocity 

most  probable  velocity  of  the  ambient  gas 


D = diameter  to  length  ratio  of  the  sample  tube 
Of  = attack  angle. 

Tabulated  values  and  the  derived  theory  for  the  R(S,  D,  Of'  function  are  given  in 
reference  3.  Ambient  pressure  is  related  to  ambient  density  by  means  of  the  ideal 
gas  law: 


1 

T“ 

a 


P 


a 


(2' 
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where 

p , p , T = ambient  density,  pressure,  and  absolute  temperature, 
a a a respectively 

M = mean  molecular  mass  of  ambient  gas 

r = gas  constant, 
m 

Although  ambient  temperature  appears  in  Eq.  (2).  the  use  of  a model  compu- 
tation for  this  parameter  is  not  as  critical  as  would  first  appear.  Subst.tuting 
Eq.  (11  into  Eq.  (2)  gives: 

M 1 d / ^a  1 (3) 

pa  = -r;  t;  pg  r<s.  n.  «r  • 

With  D and  « fixed,  R is  a linear  function  of  S at  large  speed  ratios,  as  occur  in 
satellite  night.  Thus, 

(41 

R = a S + b 


1 2 R T 
m a 


where  a and  b are  constants  and  the  other  symbols  remain  as  previously  defined. 
Hence,  the  dependence  of  the  density  calculation  upon  ambient  temperature  is  less 
than  a' square  root  dependence;  knowledge  of  Ta  is  still  needed  to  determine  s and, 
consequently,  the  sampling  function.  However,  that  dependence  can  be  made  less 
sensitive  to  changes  in  Ta  with  the  proper  choice  of  sampling  geometry,  which  will 

determine  the  coefficient  of  s in  Eq.  (41. 

In  practice  Eq.  (11  is  modified  with  an  additive  term  becoming: 


Pa  R(S,  D,  a)  + Pout 


where  P t = residual  gauge  pressure  and  the  other  symbols  remain  as  defined  in 
Eq  (11  The  residual  pressure  is  the  pressure  within  the  gauge  volume  that  does 
not  vanish  as  the  gauge  rotates  through  wake  and  enters  the  ram  phase  of  the  spin 

cycle. 
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|.|  Owgen  Rpronibi  nation 

In  calculating  the  internal  gauge  pressure,  we  have  assumed  atomic  oxygen 
recombining  to  form  molecular  oxygen.  This  is  considered  a reasonable  assump- 
tion in  view  of  the  experimental  results  obtained  from  mass  spectrometer  measure- 
ments. 4 * 6 7"8 9  The  influence  of  recombination  has  been  treated  in  some  detail  bv 
Repnev,"  and  his  results  agree  with  the  approximation  given  below.  We  assume 
that  atomic  oxygen  entering  the  gauge  recombines  to  form  molecular  oxygen  before 
ionization.  Equating  the  number  of  particles  entering  and  leaving  the  gauge  results 
in  the  following  expression: 


n (o)  v (o)  R(S,  D,  »)  A 
a a 


Vr 


n (o„)  v (o„)  A 

g 2 g 2 

2-JtT 


where 


n (o)  = 
a 

n (o_) = 
3.  2 

v (o),  v (o„)  = 
a g l 


R(S,  D,  a) 
A 


number  density  of  ambient  atomic  oxygen 
number  density  of  molecular  oxygen  in  the  gauge 

most  probable  velocity  of  ambient  atomic  oxygen 
and  gauge  molecular  oxygen,  respectively, 
tubular  factor 
tube  aperture  area. 


Solving  for  n 


g: 


n 


S 


D.  <v>) 


(7) 


(8) 
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and  2 M(o) 


M(o.^),  then 


n = 


g 


D,  a) 


Since  p = nkT,  Eq.  <6)  in  terms  of  pressure  is: 


P 


g 


R(S,  D,  a) 


(9) 


(10) 


Thus,  for  the  case  of  total  recombination,  the  internal  gauge  pressure  is  increased 
by  a factor  equal  to  the  square  root  of  two.  Introducing  a factor  for  recombination 
effects  into  Eq.  (6)  gives  the  following  expression: 

P = 
g 

where 

6 = recombination  factor;  for  no  recombination  6 equals  unity  and 
for  total  recombinations  equals  vr. 

The  above  expression  is  used  to  reduce  the  data  in  the  manner  discussed  in  the 
following  section. 


R 


Pa  R(S,  D,  «)  6 + PQut 


(11) 


5.  DATA  REDUCTION 

The  first  step  in  processing  the  ionization  gauge  measurements  was  to  obtain 
internal  gauge  pressure  measurements  as  a function  of  attack  angle  for  each  spin 
cycle.  Following  this  a linear  least  squares  fit  of  the  form 

P = m R(S,  D,  a)  + c (12) 

g 

was  performed  twice  per  spin  cycle,  once  for  values  entering  ram  and  again  for 
values  leaving  ram.  An  example  of  a fit  obtained  during  revolution  59  is  shown  in 
Figure  6.  The  values  of  the  parameters  used  in  Eq.  (9)  for  both  the  into  and  out  of 
ram  fits  are  restricted  to  attack  angles  in  the  interval  from  70  deg  to  20  deg.  This 
interval  is  further  reduced  at  low  altitudes,  due  to  the  occurrence  of  gauge  satura- 
tion. In  that  event,  the  fit  interval  is  compressed  in  the  direction  away  from  ram 
to  whatever  value  is  required  to  exclude  the  saturation  level.  The  slope  of  the 
least  squares  fit  involving  P^  and  R is: 
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Figure  6.  I.east-Squares-Fit  Data 
Acquired  on  Orbit  59.  The  166.  7 km 
altitude  is  the  altitude  at  which 
maximum  ram  occurred.  The  solid 
line  is  a least-squares  fit  to  the  data 


(13) 


From  this  equation  two  values  of  ambient  pressure  are  computed,  one  from  data 
obtained  as  the  gauge  approaches  ram  and  the  second  as  it  moves  away  from  ram. 
For  each  spin  cycle,  a mean  pressure  value  defined  in  Eq.  (14)  is  also  calculated: 


AVG 


(P  ). 
a in 


(P  ) . 

a out 


i/- 


(14) 


where  (P  ) in  and  (P  ) out  are  the  ambient  pressure  associated  with  the  into  and 
out  of  rams  slopes  respectively.  The  altitude  assigned  to  the  mean  pressure  found 
above  is  the  altitude  of  maximum  ram  (that  is,  minimum  attack  angle).  In  order 
for  Eq.  (10)  to  be  valid,  the  satellite  spin  period  must  be  sufficiently  short  such 
that  gauge  temperature,  ambient  temperature,  ambient  pressure,  and  the  recom- 
bination factor  will  not  vary  significantly  over  the  fit  interval. 

The  variation  was  tested  by  generating  the  appropriate  Jacchia  1971  (J71) 
model  values  over  spin  cycles  at  altitudes  of  300  km,  200  km,  and  perigee.  The 
variation  of  ambient  pressure  over  a quarter  spin  cycle  (that  is,  20  deg  to  70  deg 
attack  angle  interval)  for  the  model  calculated  values  was  less  than  2.  5 percent  at 
300  km,  1,  8 percent  at  200  km,  and  0,  10  percent  at  perigee.  Two  model  values 
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computed  under  the  conditions  appropriate  to  the  45  deg  attack  angle  occurrence 
entering  and  leaving  ram  were  averaged  and  the  result  compared  with  the  model 
value  computed  directly  for  maximum  ram  occurrence.  The  agreement  between 
the  two  results  was  within  0.  10  percent  over  the  altitude  region  from  300  km  to 
perigee.  Ambient  temperature  variation  over  a spin  cycle  as  determined  from 
model  computations  was  also  less  than  0.  10  percent  over  the  same  altitude  region, 
300  km  to  perigee.  The  small  variations  in  the  above  parameters  are  not  surpris- 
ing, as  for  this  flight  the  time  span  for  an  attack  angle  variation  of  50  deg  was  less 
than  2 sec. 

Before  ambient  pressures  could  be  computed  from  Eq.  (10)  a method  had  to  be 
developed  for  assessing  the  value  of  delta,  the  recombination  factor.  This  was 
performed  in  the  following  manner:  For  each  slope  determination  an  adjustment 
was  performed  wherein  delta  (6)  was  taken  as  equal  to  the  square  root  of  the  ratio 
of  the  molecular  weight  of  oxygen  to  that  of  the  mean  molecular  weight  of  the  am- 
bient gas  as  given  in  the  J71  model  atmosphere.  10,  11  At  higher  altitudes  (above 
300  km)  a lower  limit  of  16  would  have  to  be  imposed  upon  the  .171  mean  molecular 
weight.  This  adjustment,  of  course,  assumes  all  atomic  gas  species  recombine 
inside  the  gauge  and  that  the  internal  gas  is  predominantly  molecular  oxygen.  At 
perigee  altitudes  where  the  nitrogen  content  of  the  atmosphere  is  significant,  the 
numerator  of  the  delta  ratio  should  be  computed  to  include  the  molecular  weight  of 
nitrogen  as  well  as  oxygen;  however,  this  was  not  undertaken  in  the  calculation  of 
delta.  At  an  altitude  of  160  km  the  effect  of  including  molecular  nitrogen  into  the 
delta  computation  was  to  reduce  its  value  by  less  than  3 percent.  For  the  purpose 
of  this  report,  a 3 percent  change  was  not  considered  large  enough  to  justify  the 
additional  analysis  procedures  that  would  be  required  for  an  inclusion  of  a model 
composition  factor.  Also,  it  was  desirable  to  formulate  a delta  factor  that  is 
amenable  to  later  modification,  particularly  as  current  atmospheric  models  do  not 
reflect  the  results  obtained  from  more  recent  mass  spectrometer  measurements. 

5. 1 Processing 

It  was  stated  earlier,  in  the  abstract,  that  the  reduction  technique  simplified 
and  reduced  large  data  files  to  a more  manageable  size  and  format.  To  illustrate 
this,  consider  a perigee  pass  of  20  min  duration.  At  a spin  rate  of  5 rpm  a total 
of  100  spin  cycles  would  be  accumulated.  Processing  only  ram  angles  (between 

10.  Jacchia,  L.G.,  Slowey,  J.  W.  , and  von  Zahn,  U.  (1975)  Latitudinal  Changes  of 

Composition  in  the  Disturbed  Thermosphere  From  ESRO  4 Measurements, 
Smitnson.  Astrophys.  Obs.  Report  No.  309. 

11.  Jacchia,  L.G.  (19/1)  Revised  Static  Models  of  the  Thermosphere  and  Exosphere 

With  Empirical  Temperature  Profiles,  Smithson.  Astrophys.  Gbs.  fl°Dort 
No7  3IT2.  
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± 90  deg  of  the  velocity  vector)  would  produce  a file  of  9000  measurements.  One 


computer  pass  listing  of  these  point  by  point  measurements  would  require  128  pages. 


The  technique  described  in  this  re  *t  reduces  that  pass  listing  size  from  128  pages 
to  a total  of  2 pages,  without  a sig  ant  loss  of  data  quality  or  accuracy.  As  a 
result  a data  base  is  created  that  is  easily  manageable  and  adaptable  to  further 
computer  analysis. 

The  unified  data  processing  system  (DPS)  that  was  developed  is  illustrated  in 
Figure  7.  In  this  figure  the  data  flow  is  depicted  from  its  origin  onboard  the 
spacecraft  to  its  termination  point,  the  creation  of  a data  base.  The  satellite  data 
transmissions  are  pulse  code  modulated  (PCM)  and  telemetered  to  tracking  stations 
where  the  signals  are  recorded  on  tape.  The  station  raw  tapes  are  th“n  forwarded 
to  a central  processing  station  vhere  they  are  processed  and  put  into  a computer- 
compatible  format.  Referring  to  Figure  7,  the  data  flow  appears  in  four  distinct 
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Figure  7.  Functional  Flow  Diagram 


branches.  The  first  branch  is  a routine  that  processes  the  spacecraft  attitude 
information,  performs  time  correlations,  and  relates  spacecraft  attitude  to  pay- 
load  sensor  position.  The  second  branch  examines  the  experimental  data  outputs, 
unpacks,  sorts,  quality  checks  and  edits  experimental  outputs,  and  converts  and 
correlates  these  outputs  to  universal  time.  The  third  branch  is  a routine  that 
receives  tracking  data  as  an  input  and  generates  an  ephemeris  file  containing  all 
pertinent  position  parameters  as  well  as  the  appropriate  geomagnetic  parameters 
associated  with  that  position.  The  fourth  branch  is  an  indices  file  created  for  use 
in  the  computation  of  atmospheric  model  (J71)  values.  This  file  contains  quantities 
such  as  solar  declination,  K^,  ^ indices,  and  other  geophysical  parameters. 

The  density  gauge  processing  routine  requires  the  usage  of  the  data  bases 
created  in  the  four  branches  as  inputs.  Figure  8 is  a functional  flow  diagram  of 
the  density  gauge  processor.  All  necessary  initialization  occurs  in  the  first  phase 
of  the  program.  Frames  of  raw  data  are  then  input,  and  instrument  parameters 
are  calculated  on  a point  for  point  basis.  The  attitude  file  is  interrogated  and 
instrument  attack  angles  are  computed.  The  data  are  then  separated  into  the  de- 
creasing and  increasing  attack  angle  portions  of  the  ram  cycle.  The  instrument 
sampling  function,  R(S,  D,  a),  is  computed  and  gauge  pressure  is  fit  as  a function 
of  R(S,  D,  a)  for  each  spacecraft  spin  cycle.  The  slopes  of  the  fitted  curves  are 
determined  and  a value  of  ambient  pressure  is  obtained  for  both  the  into  and  out  of 
ram  portions  of  the  spin  cycle.  These  two  pressures  are  then  averaged.  Ephemer- 
is parameters  are  then  computed  and,  finally,  atmospheric  neutral  density  and 
model  density  are  calculated.  After  all  raw  data  has  been  processed,  the  density 
gauge  data  base  is  created  and  the  associated  listing  and  plots  are  generated. 


6.  GENERAL  OBSERVATIONS 

There  are  several  general  observations  worth  mentioning  in  regard  to  the  data 
processing  technique  described  in  this  report. 

(A)  The  tube  sample  function  R(S,  D,  a)  is  determined  with  the  assumption  of 
diffuse  reflection  of  particles  entering  the  sampling  tube.  If,  however,  the  par- 
ticles undergo  specular  reflection,  then  the  R factor  would  increase  in  proportion 
to  the  amount  of  specular  reflection.  At  large  speed  ratios  and  small  attack  angles 
some  specular  reflection  is  possible,  and  measurements  under  these  conditions 
would  result  in  lower  density  values  being  obtained.  The  slope  analysis  technique 
avoids  this  difficulty  by  restricting  the  fit  intervals  to  angles  greater  than  20  deg 
of  maximum  ram.  However,  although  not  shown  here,  results  that  were  processed 
at  small  attack  angles  did  not  produce  conclusive  evidence  of  specular  reflection. 


20 


PROGRAM 

INITIALIZATION 


INPUT  RAW 


TLM  DATA 


DETERMINE 

RANGE 


'CALCULATE  CURRENT 
. AND  PRESSURE  . 


CALCULATE 
ATTACK  ANGLE 


EXTRACT  INTO  AND 
OUT  OF  RAM  DATA 


CALCULATE 


DETERMINE 


CALCULATE 


CALCULATE 

EPHEMERIS 


CALCULATE 

RHO 


CALCULATE 


MODEL 


IDG 

DATA 

BASE 


PLOT 


Figure  8.  Systems  Flow  Diagram 


(B)  As  instrument  response  is  processed  as  a function  of  time,  undetected 
variations  in  spin  axis  position  will  have  their  greatest  impact  at  small  attack  angles 
as  the  spin  period  is  unaffected  by  the  direction  of  the  spin  axis.  For  example. 


consider  the  instrument  attack  angles  of  10  deg  and  40  deg  occurring  at  times  tj 
and  t2>  respectively;  then,  a 10  deg  spin  axis  misalignment  produces  a 40  percent 
change  in  the  attack  angle  occurring  at  time  tj,  changing  from  10  deg  to  14.  1 deg, 
whereas  at  time  the  change  is  about  3 percent  - from  40  to  41.  2 deg.  Of  course, 
this  result  is  dependent  upon  the  mode  of  orientation  and  may  not  apply  in  the  event 
of  a different  satellite  orientation. 

(C)  Another  consideration  along  with  the  two  already  mentioned  is  the  relation 
ship  of  the  R factor  to  the  gauge  sampling  geometry.  The  data  shown  in  Table  1 
were  obtained  from  reference  3 and  shows  the  sampling  function  (R(S,  D,  a))  for 
two  tubulation  ratios  as  a function  of  attack  angle.  These  sampling  functions  are 
applicable  to  cylindrical  tubes  defined  in  terms  of  their  diameter  to  length  ratio. 
However,  there  are  cases  where  a gauge  tubulation  ratio  cannot  be  precisely  de- 
fined in  terms  of  the  sampling  tube  dimensions,  resulting  in  an  equivalent  or 
effective  ratio  in  place  of  a geometric  one.  The  uncertainty  between  the  two  tubula 
tion  ratios  would  result  in  a significant  effect  at  small  attack  angles,  but  the  effect 
would  rapidly  diminish  as  the  attack  angle  increases.  This  can  be  seen  from 
Table  1 where  the  differences  in  R factors  at  0 deg  is  26  percent,  while  at  30  deg 
the  difference  is  3 percent,  and  these  differences  are  for  gauge  tubulation  ratios 
that  differ  by  40  percent. 


Table  1.  Tubulated  Gauge  Sampling  Functions 


speed  ratio,  u/v 
r m 

diameter /length  ratio  of  sampling  tube 


R(S,  Dt,  a) 

R(S,  D,  a) 

— 

Attack  Angle  (a) 

(D  = 0. 30) 

(D  = 0.50) 

(Deg) 

110.  3 

87. 7 

0 

82.  7 

74.  1 

10 

50.  3 

52.4 

20 

34.  5 

35.7 

30 

25.  0 

25.  6 

40 

18.  1 

18.  3 

50 

12.4 

12.4 

60 

7.  53 

7.49 

70 

3.43 

3.  38 

80 

0.  515 

0.  504 

90 

(D)  A final  observation  concerning  the  analysis  technique  is  its  tendency  to 
minimize  the  effects  of  certain  types  of  systematic  errors.  An  example  of  this 
would  be  a change  in  the  spin  frequency  that  produced  a mislocation  of  the  mini- 
mum ram  angle.  If  the  angular  offset  were  in  the  direction  leaving  ram,  then  the 
out  of  ram  gauge  values  would  be  processed  with  R factors  that  were  too  large; 
consequently,  a smaller  value  of  slope  and  ambient  pressure  would  be  obtained. 
Conversely,  the  R factor  values  used  to  process  the  gauge  pressures  entering  ram 
would  be  too  small.  The  subsequent  averaging  of  the  two  pressures  would  have  the 
effect  of  compensating  for  the  aspect  error.  The  same  may  be  said  of  hysteresis 
effects  present  in  the  gauge  response,  where  slight  differences  in  response  can 
occur  depending  upon  whether  the  internal  gauge  pressure  is  increasing  or  de- 
creasing. As  approximately  the  same  pressure  region  is  being  covered  during 
the  entering  and  leaving  phases  of  the  spin  cycle,  the  hysteresis  effects  will  also  be 
reduced  by  the  averaging  process. 


T.  DATA  RESULTS 

Data  shown  in  Appendices  A and  B were  obtained  during  the  period  from  5 
November  1974  to  17  January  1975.  The  250  orbits  shown  are  not  all  inclusive; 
additional  orbits  from  this  period  are  awaiting  processing,  and  others  require 
reworking  to  remove  excessive  noise.  The  orbits  presented  here,  however,  are 
a representative  sample  of  the  data  obtained.  Shown  in  Appendices  A and  B 
are  two  sets  of  density  plots  and  a chart  of  satellite  positional  data.  Figure  A 1 
shows  the  satellite  position  in  terms  of  perigee  location  where  the  local  time, 
latitude,  and  altitude  are  given  as  a function  of  time  and  revolution  number.  In 
Appendix  A,  plots  of  density  (gm/cc)  vs  universal  time  are  shown.  The  circles  are 
values  of  density  obtained  from  gauge  measurements,  and  the  crosses  are  the 
appropriate  Jacchia  71  model  values  (calculated  under  the  same  conditions  as  the 
measurements).  In  the  upper  portion  of  the  same  plots,  the  crosses  denote  the 
ratios  of  the  measured  densities  to  the  J71  model  values.  The  shade- sun  nomen- 
clature shown  in  the  legend  represents  the  time  (U.  T.  ) when  the  satellite  was  in 
darkness  and  sunlight,  respectively.  Other  designations  shown  in  the  ordinates, 
abscissa,  and  legend  are  self-explanatory. 

The  ratios  of  measured /model  density  values  are  generally  higher  than  unity. 
This  is  partially  due  to  the  method  of  processing  gauge  pressures  wherein  the  sen- 
sitivity factor  for  nitrogen  gas  was  used  and,  therefore,  the  results  are  given  in 
terms  of  equivalent  nitrogen  pressures  that  have  subsequently  been  converted  to 
mass  densities.  As  the  atmosphere  consists  of  a mixture  of  gases,  the  sensitivity 
factors  for  the  individual  gas  species  should  properly  be  included  in  the  pressure 


conversion  process.  However,  to  accomplish  this  the  relative  concentrations  of 
each  gas  specie  must  be  determined.  Hence,  for  the  purposes  of  this  report, 
such  processing  was  not  attempted  as  it  was  felt  that  the  use  of  equivalent  nitrogen 
pressures  provided  a more  consistent  data  base  that  can  later  be  modified  to 
accommodate  compositional  data.  The  relative  sensitivity  factor  for  molecular 
oxygen  (R  = Sq^/Sj^)  is  1.  09.  Thus,  if  the  atmosphere  were  total  oxygen,  the 
pressures  would  be  correspondingly  reduced  by  9 percent.  Gauge  sensitivity 
factors  and  their  application  will  be  discussed  in  more  detail  in  a second  report 
covering  the  data  obtained  during  the  period  from  February  to  May  1975. 

In  the  altitude  plots  given  in  Appendix  B,  the  circle  and  cross  symbols  denote 
gauge  measurements,  where  the  circles  are  for  measurements  performed  when 
the  satellite  descended  in  altitude  towards  perigee  and  the  crosses  are  for  measure- 
ments performed  when  the  satellite  ascended  in  altitude  away  from  perigee.  The 
solid  lines  shown  are  least  squares  fits  through  the  data  points.  In  some  altitude 
plots  further  processing  is  required  to  remove  noise  values.  It  is  readily  seen 
from  examination  of  the  density  vs  time  plots  which  values  are  noise,  and  refer- 
ence to  the  corresponding  altitude  plot  will  show  the  distortion  the  noise  value  has 
upon  the  least  squares  fit.  At  higher  altitudes  the  downleg  (circles)  and  upleg 
(crosses)  profiles  diverge.  A major  portion  of  this  divergence  is  due  to  the  spatial 
separation  of  the  measurements.  Perigee  is  initially  located  at  midlatitudes  in  the 
northern  hemisphere  and  precessing  northward  to  the  pole.  The  downleg  measure- 
ments are  performed  at  high  northern  latitudes,  whereas  the  upleg  measurements 
are  taken  at  middle  and  equatorial  latitudes  nearer  to  the  sub- solar  point.  The 
upleg  measurements  therefore  reflect  the  higher  densities  associated  with  the  larger 
temperatures  due  to  the  proximity  of  the  solar  bulge. 
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